Impurity and interface bound states in d x 2_ y 2 + id xy and p x + ip y superconductors 
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Motivated by recent discoveries of novel superconductors such as Na^CoC^ • 2/H2O and S^RuCU, 
we analysize features of quasi-particle scattering due to impurities and interfaces for possible gapful 
d x 2^ y 2 + id xy and p x + ip y Cooper pairing. A bound state appears near a local impurity, and a 
band of bound states form near an interface. We obtained analytically the bound state energy, 
and calculated the space and energy dependent local density of states resolvable by high-resolution 
scanning tunnelling microscopy. For comparison we also sketch results of impurity and surface states 
if the pairing is nodal p- or d-wave. 



Recently, Takada et al discovered a novel supercon- 
ductor Naj;Co02 • J/H2O (x = 0.35) with a supercon- 
ducting transition temperature T c = 5K. [1] A few fea- 
tures of this material bare strong connection to cuprates: 
1) It has a layered structure. 2) As Cu 2+ in cuprates, 
Co 4+ atom is in a spin- 1/2 state according to first prin- 
ciple calculation of Singh. [2] Combined with the fact 
that the cobalt triangular lattice is frustrating to anti- 
ferromagnetic ordering, the new material offers a likely 
situation for the physics of Anderson's resonating valence 
bond (RVB) theory. [3] Soon after the discovery, theories 
based on RVB physics [4-7] and renormalization group 
analysis [8] predicted d+id'-w&ve pairing (d = d x 2_ y i and 
d' = d xy ), while other theories suggest p x + ip y -w&ve pair- 
ing derived from the weak ferromagnetic instability, [9] in 
close analogy to the case of Sr2RuO4.[10] Identifying the 
pairing symmetry would be a necessary step toward the 
understanding of the new superconductor. In this paper, 
we propose tunnelling measurements of impurity and in- 
terface states that are sensitive to both the gap amplitude 
and the internal phase of the gap function. Such measure- 
ments have played invaluable roles in high temperature 
superconductors in the context of nodal d-wave pairing 
.[11, 12] Our main results are as follows. As a conse- 
quence of the full gap as well as the internal phase degrees 
of freedom of d + id' and p x + ip y Cooper pairs, a bound 
state appears at any nonzero scattering intensity near a 
local impurity and a band of bound states form near an 
interface. The bound state energy is near the gap edge at 
weak scattering strength, and it approaches zero energy 
(the Fermi level) at increasing scattering strengths. We 
also calculated the energy and space dependent LDOS, 
whose rich features are directly resolvable by future STM 
and can help identify the pairing symmetry in the new su- 
perconductors. For comparison, we also mention briefly 
the results of impurity and interface states for nodal p- 
and d-wave pairing. 

As usual the elastic scattering problem is best de- 
scribed in terms of the retarded T-matrix formulation, 



where a, b denotes the position of the impurities and all 
other notations are standard. We suppressed the energy 
dependence in the Green's functions, as it is conserved 
in elastic scattering. The T-matrix is given by 



(T- 1 )(a,6) = (y- 1 )(a,6)-G (a,6), 



(2) 



where V(a, b) is the general impurity potential that may 
be off-diagonal. In our case, G, Go and T are further 2x2 
matrices in the particle-hole Nambu space. The scatter- 
ing problem is solved once Go is known. The LDOS at 
site i is given by 

N(i,u) = -Im[Gn(i,i;w) + Ga2(t,i;-w)]/7r, (3) 

with the energy argument ui restored. A peak in N(i, ui) 
appears if either Gn(i,i;w) or G22(m;~ w ) diverges. 
This bound/resonance state occurs if Det[T _1 (±w)] = 0. 
Due to the mixing of particle and hole in the presence of 
pairing, it is possible that there are two peaks in N(i,io) 
but Det(T _1 ) = is satisfied at only one energy, or vice 
versa. In the following discussion, we always count the 
bound/resonance states according to the peaks seen in 
the total density of states N(i,u>). 

Let us write the gap function as , in the momentum 
space, Ak = Ae J ' 6 ' k where A is the gap amplitude, 6u 
is the azimuthal angle of the vector k and I = 0, ±1, ±2 
for gapful s-, p- and d-wave pairing, respectively. We 
include the case of s-wave pairing for comparison. The 
above pairing function is of simplified form, suitable near 
the normal state Fermi surface, and suffices for qualita- 
tive discussion of low energy quasi-particle states. Then 
G (i,j) = G (r) (with r = r t - r,) is given by 
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G(i,j) = G Q ^j)+Y,Go(*,a)T(a,b)G (b,j), (1) 

a,b 



where lu + = lu + i0 + , v = ±, is the normal state 
energy dispersion, ctq is the 2x2 unit matrix, 01, 2, 3 
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are the Pauli matrices, a± = (<ti ± io~i)j2 and Ji(u) — 
/ * d9cos(W) exp(iu cos 9) /2tt is the Bessel function. In 
arriving at the above results, we have assumed a cylin- 
drical Fermi surface with Fermi vectors of magnitude hp, 
and constant density of states -/V near the Fermi level. 
We emphasize that a particle-hole asymmetry is present 
in the normal state DOS of Na^CoC^- We shall comment 
on such effects without going into details in the following 
qualitative discussions. 

I. Scattering from a local impurity: In this case 
we set the impurity site at the origin, i.e., a = b = 0, 
and drop these indices in V = V m <7 + V S <J3 and T^ 1 = 
V^ 1 — Go(0,0), where V m ,s is the strength of (classical) 
magnetic/scalar potential. With Go in Eq.(4) at hand, 
the T-matrix is now given by 

T " 1= ^ 1+ VA ^ +)2 (^ + A ^)- < 5 > 

One sees that Im(T _1 ) — > in the sub-gap regime 
lo 2 < A 2 , so that a true bound state could be generated in 
this regime provided Dct(T _1 ) = 0. This should be con- 
trasted to the case of virtually bound impurity state, or 
the resonant impurity state, in the case of nodal d-wave 
pairing. [13] The condition Dct(T _1 ) = is governed 
by the value of the dimensionless scattering strengths 
c m .s = ^N{)V m ^ s . A few cases are classified as follows. 

The case of s-wave pairing (I = 0) has been discussed 
in the literature, [14, 15] and we list briefly some of the 
results for comparison: 1) No sub-gap bound states can 
be generated for a scalar impurity (V s ^ and V m =0). 
This is consistent with the Anderson theorem that s-wave 
pairing is robust against scalar impurities. 2) A magnetic 
impurity (V s = and V m ^ 0) can always generate a 
pair of bound states at energies cob = ±|A|(1 — c 2 „)/(l + 
c 2 m ). In any case, the spatial modulation of the LDOS 
in the presence of the impurity is given by N(i,u>) — 
Af(w) oc Jg(k F r l ) from Eqs.(l), (4) and (5) under the 
given approximation. Here Af(u>) is the site-independent 
bulk density of states in the superconducting state. 

For p- and d-wave pairing, the off-diagonal <ji- 
componcnt in T^ 1 is zero. This is not an accidental 
result from the adopted approximation, but rather a rig- 
orous result from the pairing symmetry, which forbids 
the on-site pairing amplitude (related to the anomalous 
part of Go(0)) to be finite. Consequently, both scalar 
and magnetic impurities can generate bound states. 1) 
For a scalar impurity, Dct(T _1 ) = is satisfied at 

uo b = ±|A|/vT+# (6) 

In general this implies two peaks in the LDOS according 
to Eq.(3). However, depending on the ratio to 2 / A 2 one of 
the peaks may dominate over the other, with an associ- 
ated change in the spatial dependence of the LDOS. We 
present a few examples of the energy and space depen- 
dent LDOS in Figs.l for Cooper pairing with I = 1,2. 




FIG. 1: Density of states as a function of energy lo and the 
radial distance r off a scalar impurity. See the text for details. 




FIG. 2: The same plot as Figs.l but for a magnetic impurity. 

For the weak impurity case c s = 0.5 in Figs.l (a) and 

(c) , ujb is near the gap edge, the dominant peak is at the 
energy with opposite sign to c s , and the corresponding 
DOS right at the impurity site is maximal. In contrast, 
for the strong impurity case c s = 10 in Figs.l (b) and 

(d) , LUb is approaching the Fermi level (zero energy), the 
dominant DOS peak energy has the same sign as that of 
c s , and the corresponding DOS is vanishing right at the 
impurity site. Note that the cusp-like feature at lo = ±A 
away from the impurity is just a feature of Af(uj). 2) For 
a magnetic impurity, Det(T _1 ) = is satisfied at 

u b = -sign(c m )|A|/Vl + c 2 „- (7) 

Since T _1 cx ero in this case, there are actually two peaks 
in DOS, according to Eq.(3), located symmetrically with 
respect to the Fermi level. Examples are shown in Figs. 2, 
in comparison to Figs.l. By inspection, we see that ex- 



3 



cept for the symmetrical peaks, Figs. 2 are basically sim- 
ilar to Figs.l. On the other hand, in both scalar and 
magnetic impurity cases the difference between I = 1 
and I = 2 is mild. This would pose difficulty for STM 
to resolve this quantum number. Fortunately this can 
be resolved easily by other means such as spin suscep- 
tibility measurements from the fact that singlet pairing 
(1 = 2 here) forms a gap for spin excitations while triplet 
pairing (1 = 1 here) does not. 

It is pertinent at this stage to comment on the effect of 
particle-hole asymmetry in the normal state Fermi sur- 
face. As can be seen from the derivation of Go, this would 
introduce a 03-component in Go (0,0), which effectively 
acts as an excess energy- dependent scalar potential in 
T^ 1 . Therefore, the effect is to modify the bound state 
energy, and to break the symmetry of the bound state 
energies in the case of magnetic scattering. 

II. Scattering from an interface: We shall model 
an interface by an extended line of impurities. This could 
be fabricated by chemical erosion. It is to our advantage 
in that the T-matrix formalism can still be applied. Since 
the unperturbed system at hand has rotational symme- 
try, the interface states should not depend on the surface 
normal direction n, which we fix to be h = y for def- 
initcness. Due to the remaining translation symmetry 
along the x-axis, we can do partial Fourier transforms of 
Eqs.(l), (2) and (4) with respect to x to find the reduced 
t-matrix equations at the x-direction wave vector q as 

g(yi,yj) = go(yi - Vj) + go(yi)tgo(-yj), 

2irN (u+ cospyao + A J2 U cos(py + vlO q )a v 



+ 
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-p^A 2 - uj\ 
= (oj + (Jo + A cos 16 \o\ ) , 



where v = V m ao + V s a^ is the same as the form of 
a single impurity, p = \Jk 2 F — q 2 = kp\sm9 q \ and 
6 q = arccos(<7/fcp). The conserved momentum q is sup- 
pressed in the arguments of g, go and t~ x for brevity. 
The problem is reduced to an effective single impu- 
rity scattering in one-dimension. With the implicit w 
and q arguments restored, the partial DOS is given by 
N(w,y;q) = -lm[g n (u), y; q) + g 2 2(~u), y; q)]/n, and the 
total density of states is N(u,r) = J dqN(u,y; q)/2w, 
which is independent of x due to the translation symme- 
try. (The integration over q should be cutoff at ±fcp.) 
Again Dct(f _1 ) = would predict a bound state. 

Although the above formulation is versatile to deal 
with any value of V m . s , we shall consider only the more 
likely scalar interface with V m = 0. It is easy from the 
above equations that bound states occur at energies given 

by 



IAc 2 i 



kp sin 2 q 



Ac 2 + hi sin 2 



(8) 
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FIG. 3: Dispersion of the positive bound state energy as a 
function of the wave vector along the interface. See the text 
for details. 



(a) 1=1, c =0.5 




(b) 1=1, c =10 




(c) l=2, c =0.5 



(d) l=2,c=10 




FIG. 4: Density of states as a function of energy uj and the 
distance d off an interface. See the text for details. 



which clearly form two bands. (Note that we have taken 
the lattice constant to be unity so that kp is dimension- 
less.) It is also clear that no sub-gap bound states exist 
for s-wave pairing (/ = 0). 

The dispersion of the positive bound state energy is 
plot in Fig. 3 for weak (c s = 0.5, thick lines) and strong 
(c s = 10, thin lines) interface with p-wave (1 = 1, solid 
lines) and d-wave (1 = 2, dashed lines) pairings. Here 
we have set kF — tt/2 for calculation. One sees that 
the energy disperses near the gap edge for weak scatter- 
ing interfaces, and it tends to cover the whole sub-gap 
regime for strong interface scattering. Furthermore, the 
difference between p- and d-wave pairing is reflected in 
the number of minima, being identical to I, in the disper- 
sion. 

The spatial dependence of the LDOS near the interface 
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can be calculated from the above theory. Examples are 
shown in Figs. 4. Consistent with the above bound state 
energy dispersion, the sub-gap states are near the gap 
edge (or tend to cover the whole gap regime) for a weak 
(or strong) scattering interface, distributed more or less 
symmetrically (or asymmetrically) with respect to the 
Fermi level. In the limit of unitary scattering c s — > oc 
(not shown here) the LDOS becomes symmetrical in en- 
ergy again. An interesting feature in Figs. 4(b) and (d) is 
that the peaks or bumps in energy oscillate with increas- 
ing distance from the interface, forming wave-like pattern 
in the energy-distance space. It is also clear from Figs. 4 
that the spatial profile decays much more slowly than the 
single impurity cases in Figs.l and 2. 

We note that Matsumoto and Sigrist [17] have ad- 
dressed the quasi-particle states near a sample surface 
and a topological domain wall (with a n phase shift in 
the pairing gap) for p x + ip y -w&vc pairing in terms of 
quasi-classical theory. They found that sub-gap states 
appear near the domain wall but not the surface. Our 
interface is actually a non-topological domain wall but 
with potential scattering. 

III. The case of nodal p- and d-wave pairing: Along 
similar lines to that sketched above, we have also consid- 
ered the impurity and interface states for nodal p- and 
d-wave pairing for comparison. The gap function may 
be written as Ak = AsinZ#k or Ak = AcosZ#k, depend- 
ing on whether one of the nodal or antinodal directions 
is along the x-axis. Note that there are only one nodal 
and one antinodal direction for p-wave pairing. Due to 
limited space we sketch the results without going into 
details. 

For the local impurity case, we found resonant energies 
at io r ~ ±7rA/[Zc TOiS ln(4Zc TO;S /7r)] for a scalar/magnetic 
impurity. This reduces to the known result in the case of 
nodal d-wave pairing (/ = 2). [13] The new features for the 
case of nodal p-wave pairing is that that LDOS pattern 
near the impurity is two-fold symmetric, forming stripe- 
like features extending along the anti-nodal direction, 
which should be compared to the four-fold symmetric 
pattern in the case of nodal d-wave pairing. [13, 16, 18, 19] 

On the other hand, the interface states depend on 
the interface orientation: 1) If the (scalar scatter- 
ing) interface is along one of the nodal directions, 
say x, there are bound states at energies uJb = 

± A k F sin 9 q sin W q \/ y 1 Ac 2 s + k\ sin 2 9 q . The definition of 
q is the same as in section II. These energies all approach 
zero in the unitary limit c s — > oo. The abundance of zero 
energy states is due to the fact that in this limit quasi- 
particles reflect spectacularly from the interface, experi- 
encing a sign change of the gap. The same physics is 
nicely described in Refs.[20, 21] in other contexts. 2) 
Finally if the interface is along one of the anti-nodal di- 
rections, redefined also as x, there are resonant states 
exactly at ui r = ± A cos I6 q irrespectively of the scatter- 



ing strength. In fact this is equivalent to the case of 
s-wave pairing but with a g-dependent gap amplitude. 

IV. Closing remarks: We have only shown results for 
the cases I — ±1,±2 that are relevant in the new su- 
perconductors, but the theory is clearly general for any 
integer value of I. There are some details missing in the 
theory, however. First, it does not take into account pos- 
sible anisotropy in the normal state Fermi surface. For 
example, in Na a; Co02, the Fermi surface has a rounded 
hexagonal structure. [2] Such anisotropy may cause cor- 
responding anisotropic LDOS pattern around impurities. 
Second, Eq.(4) is obtained by fixing the momentum on 
the Fermi surface while integrating over energy. This 
possibly leaves out an excess decay of Go in space with 
the length scale £ = vp/A. Apart from such details, our 
qualitative analytical results are robust. 
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